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ABSTRACT 

We have selected a sample of local E+A galaxies from the Sloan Digital Sky Survey 
(SDSS) Data Release 7 for follow up integral field spectroscopy with the Wide Field 
Spectrograph (WiFeS) on the ANU 2.3-m telescope. The sample was selected using 
the Ha line in place of the [OIIJA3727 line as the indicator of on-going star formation 
(or lack thereof). This allowed us to select a lower redshift sample of galaxies than 
available in the literature since the [OIIJA3727 falls off the blue end of the wavelength 
coverage in the SDSS for the very lowest redshift objects. This low redshift selection 
means that the galaxies have a large angular to physical scale which allows us to resolve 
the central ~1 kpc region of the galaxies; the region where stellar population gradients 
are expected. Such observations have been difficult to make using other higher redshift 
samples because even at redshifts z ~0.1 the angular to physical scale is similar to 
the resolution provided by ground based seeing. Our integral field spectroscopy has 
enabled us to make the first robust detections of Balmer line gradients in the centres 
of E+A galaxies. Six out of our sample of seven, and all the galaxies with regular 
morphologies, are observed to have compact and centrally-concentrated Balmer line 
absorption. This is evidence for compact young cores and stellar population gradients 
which are predicted from models of mergers and tidal interactions which funnel gas 
into the galaxy core. Given the generally isolated nature of our sample this argues 
for the galaxies being seen in the late stage of a merger where the progenitors have 
already coalesced. 

Key virords: galaxies: evolution - galaxies: formation - galaxies: stellar content 



1 INTRODUCTION 

E+A galaxies, their optical spectra characterized by strong 
Balmer absorption lines and little or no emission, repre- 
sent galaxies undergoing rapid evolution in their star for- 
mation properties. The strong Balmer lines are indicative of 
a substantial population of A-type stars which must have 
formed within the last -^1 Gyr. However, the lack of opti- 
cal emission lines imply that the episode of star formation 
in which they formed has ended and star formation is not 
ongoing. The normal interpretation of this spectral signa- 
ture is that the galaxy is being observed in a post starburst 
phase (e.g. ICouch fc Sharpled Il987l : iPoggianti et~aLl Il999l ) 
although it can also be produced by the abrupt trunca- 
tion of star formation in a disk, without necessarily requir- 



ing a starburst f e.g. IShiova et all |2004 iBekki et al] l2005l : 
IPracv et al.ll2005h . This spectroscopic evolution may be ac- 
companied by a morphological one where a star-forming 
disk galaxy is transformed i nto a quies cent spheroidal sys- 
tem (C aldwel l et al. 1996; Za bludoff et al. 1996). 



Much attention has been paid to the environmental 
influences which give rise to the abrupt changes in star 
formation rate within E+A galaxies and several processes 
have been suggested . These include major and unequal 
mass galaxy mergers (iMihos fc Hernquistlll996l : iBekki et al.l 
gal axy-galaxy intera ctions with mass ratios close 
to unity (|Bekki et all |2005| ) and, in the clus ter environ- 
ment, interactions with the cluster tidal field (|Bekki et al.l 
[2001j), galaxy harassment (|Moore et al., ,19981 ') and interac- 
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tion with the intra-c l uster medium |"Gunn fc GottI 1 19721 : 
IDressler fc Gunnll 19831 : iBothun fc Dressleill 19861 '). While the 
situation in gala-x y clusters remains somewhat unclear 
l|Pracv et al.ll20ld ). there is growing consensus that galaxy 
mergers are the mechanism responsible for producing the 
local field E+A population as seen in the analysis of E+A 
samples from large redshift surveys. Ground based imaging 
and HST follow up of a sample of 21 E+A galaxies selected 
from the Las Campanas Redshift Survey (LCRS) revealed an 
increased rat e of tidal features implying galaxy merg ers had 
taken place (|Zabludoff et al.lll996l : IVang et al.ll2008l i. Using 
a sample selected from the Two Degree Fi eld Galaxy Red- 
shift Survey (2dFGRS'). iBlake et al.[ (|2004 ) concluded that 
the significant incidence of tidal disruption, the nature of 
the luminosity function, and the results of a near-neighbor 
analysis, wer e cons i stent with them being in the late stage 
of a merger. iGotol (Hooi) found a higher incidence of near 
neighbors for E-l-A galaxies in the SDSS suggesting galaxy- 
galaxy interactions may also be an important mechanism. 

Another, more direct, way of discerning the physical 
cause of the E+A phase is to examine the internal proper- 
ties of the galaxy itself. The internal kinematics and spa- 
tial distribution of the young and old stellar populations 
are expected to differ between different formation scenar- 
ios. For example, in a major merger the young stellar pop- 
ulation should quickly become dynamically p ressure sup- 
ported l|Norton et al.l l200ll : iBekki et al] 120051 ). For more 
extreme mass ratio mergers and tidal interactions, signif- 
icant rotation of the young st ellar population is gener- 
ally expected (|Bekki et all 120051 ). The two major spatially 
resolved spectroscopic studies of low redshift 'field' E+A 
galaxies (|Norton et al.|[200ll : |Pracv et al.ll2009l ) have reached 
quite different conclusions on the prevalence of rotation in 
E+ As. Using follow-up slit spectroscopy of the LCRS sam- 
ple, iNorton et al.l (|200l|) f ound Uttle evidenc e for rotation 
in most cases. In contrast. iPracv et al.l l|2009[ l's two dimen- 
sional Integral Field Unit (IFU) spectroscopy of a sub- 
sample of the iBlake et al.) \200i ] 2dFGRS E+A galaxies 
showed rotation of the young stellar population to be ubiq- 
uitous. This discrepancy could partly be explained by the 
rotation axis not always corresponding to the isophotal ma- 
jor axis, dampening th e amount of rotati on measured using 
a one dimensional slit (|Pracv et al.ll200^ ). 

If the E+A phase is the result of gas-rich mergers, 
then evidence of this should be contained in the spatial 
distribution of the stellar populations. During a merger, 
tidal torq ues transfer angular mome ntum from the gas to 
the stars ([Barnes fc Hernguistl [l99^ ). driving the gas to- 
ward the centre and result i ng in a compact central starburst 
dMihos fc Hernguistl [1993 . 1 19961 : [b arnes fc Hernguistl 1 19961 : 
iHopkins et al."2009'). Feedback processes then halt the star- 
burst (Springcl & Hcrnguist 2003: Sp ringel et al . 2005 ) leav- 
ing an old stellar populatio n that is distributed like a nor- 
mal e arly-type galaxy (e.g. lHernguis3ll992l : IHopkins et al.l 
I2OO9D and a you ng component that forms a centrally- 
concontrated cusp (|Mihos fc Hernguist|[l993 : iHopk ins et al.l 
L2OO9). This young c omponent should be quite c o mpact , with 
scale s of <1 kpc (iBarnes fc Hernguistl 1 19961 : iBekki et al.l 
I2OO5I : IHopkins et aUbOOgi r and is a possible explanation of 
the excess light observed in the cent r es of merging star- 
burst galaxies (iHibbard fc Yunl 1 19991 : iRothberg fc Josephi 
l2004 ). iKohno et all (|20oi r detected centrally concentrated 



diffuse and gravitationally stable molecular gas in the post- 
starburst galaxy NGC 5195 which may be the remaining 
component of the gas used in fueling the central starburst. 
In the E+A galaxies the remnant of such a starburst should 
be observable as a Balmer line absorp tion enhancement 
and gradient with in the central ~1 kpc (|Pracv et al ] |2005l : 
iBekkiet al.]|2005h . These predictions of mergers resulting 
in a centrally concentrated post-starburst region and strong 
radial gradients have been strengthened by a recent set of 
detailed fully three-dimensional hydrodynamic simulations 
of E+A galaxy production from galaxy merging taking into 
account a wide range of initial parameters such as mass, gas 
fract ion, orbital configu ration, mass ratio and dust proper- 
ties (jSnvder et aLll201ll ). 

Attempts to measure the spatial distribution of the 
post-starburst stellar population in E+A galaxies has been 
hindered by the small expected angular scale of the star- 
burst remnant at the redshi fts of E +A galaxies selected 
from large redshift surveys (B lake et al. 2004: £oto 200'^), 
where it is smaller or comparable to the typical seeing 
scale. For example, at the median redshift of these surveys 
{z ~0.1) the 1 kpc expected scale of the starburst projects 
to ~0.5 arcseconds. iNorton et al.l (|200ll ) compared the ra- 
dial distribution of the young and old populations in the 
LCRS sample and concluded that the young stellar popula- 
tion was centrally c oncentrated. The seeing FWHM of the 
INorton et al.l ([2003) observations corresponded to ~2-3 kpc 
and in about half the cases they were able to spatially resolve 
the young stellar population by showing it to be broader 
than the stellar PSF. The implication, in these cases at least 
was that, the starburst population is not confined to the 
galaxy core (~1 kpc). Likewise, follow u p of a Hl-detected 
post starburst galaxy from the SDSS bv IChilingarian et al.l 
( 2009. ) found the very young stellar p opula tion spread over 
the central ~2.5 kpc. lYagi fc Goto! (|2006f ) analyzed three 
early type E+A galaxies selected from the SDSS using long- 
slit spectroscopy. In the one galaxy they could successfully 
estimate a H5 equivalent width profile, they found that the 
post-starburst signature was as extended as the continuum 
light, and in the other two galaxies concluded (on statisti- 
cal grounds) that the po st-starburst region was also signif- 
icantly extended. iPracv et al. (2009l) were unable to detect 
any Balmer line gradients in their sample although their 
limited spatial coverage coupled with spatial resolution con- 
straints imposed by the delivered image quality meant they 
could not robustly rule out compact gradients in the stellar 
populations. 

In this paper we circumvent these physical scale resolu- 
tion issues by investigating a new, very low redshift sample 
of E+A galaxies which have large angular to physical size 
scales. This is achieved by modifying the standard definition 
of an E+A galaxy, using the Ha line as the constraint on 
emission (and hence ongoing star formation) instead of the 
traditional [OIIJA3727 line. This change allows selection of 
a local sample of E+As at 2; < 0.01 to easily resolve their 
central 1 kpc. 

In what follows. Section 2 describes the details of the 
sample selection; Section 3 describes the observations, data 
reduction and analysis; in Section 4 we present our results, 
including spatially resolved Balmer line maps; and in Sec- 
tion 5 we discuss these results. Our conclusions are given in 
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Section 6. Throughout, we adopt an Q,m = 0.3, Q,a = 0.7 
and Ho = 70knis~^ Mpc~^ cosmology. 

2 THE SAMPLE 
2.1 Selection criteria 

E+A galaxies are rare i i i the local u niverse (jZabludoff et al.l 
1 19961 : Isiake et al.]|2004l : iGotd 120071 ') and there are very few 
nearby j z < 0.01) example s known. For example, from the 
2dFGRS [Brake et al.l (|2004l ') cataloged 56 E+A galaxies but 
only two of these had redshifts of z < 0.01. The catalogue 
of E+A gal axies produced from the much larger SDSS by 
iGotd (|2007l ') has no E+As with redshifts oi z < 0.03. The 
dearth of E+As is the result of the spectral range of the 
SDSS spectra (A >3800A) not containing the [OII]A3727 
line for the very lowest redshift objects and not a result of 
a true lack of local E+A galaxies in the SDSS. Although 
the E+A fraction does increase rapidly with redshift (e.g. 
iPoggianti et"allll999l : ISnvder et aLll201ll ) . The use of low or 
non-existent equivalent widths of the [OII]A3727 line as the 
constraint on ongoing star formation is historical, since E+A 
galaxies were original ly discovered and studie d in intermedi- 
ate redshift clusters (jDressler fc Gunnlll983l ) where the Hq 
line was not available. In fact, at low redshift the Hq line is 
a more natural choice for this since it is more easily relatable 
to the star formation rate, intrinsically stronger, less prone 
to dust obscuration and has been shown to be a more robust 
const raint on star formation in E+A galaxies (|Blake et al.l 
|2004 ). By replacing the [OII]A3727 constraint with one de- 
fined using the Ha line we have selected a. z < 0.01 E+A 
sample from the SDSS data release 7 spectral line catalogue. 
The selection criteria used were (our definition of equivalent 
width has emission denoted as negative values): 

• Little or no Ha emission: EW(Ha)> -2.5 A 

• Strong absorption in the U5 line: EW(H5)> 5.0 A 

• We also demand strong absorption in the H7 and H/3 lines 
since th is has been shown to select a much higher fidelity 
sample (|Blake et al.ll2004l ): EW(H7) & EW(H/3)> 4.0 A 

• Small errors on the equivalent width measurements of the 
Balmer absorption lines: The errors in each of EW(H(5), 
EW(H7), EW(H/3)< 1.0 A 

• High quality redshift: probability the redshift is correct 

2^coiifidcncc >0.95 

• Low redshift: 0.0005 < z < 0.01 

• Sufficiently bright for IFU follow up: r < 15.5 mag 

In addition to these criteria there were also constraints on 
the right ascension and declination for follow up observa- 
tions. In total this gave us eight possible targets that could 
be observed during our scheduled telescope time and we ob- 
tained data for all but one. In the absence of right ascen- 
sion and declination constraints, twenty galaxies in the SDSS 
meet these criteria. Imposing constraints on the H/3 and H7 
lines in addition to H(5 has a significant effect on the sample 
size, particularly H/3 which is prone to emission line filling. 
Removing the H/3 constraint increases the sample size from 
twenty to thirty whilst removing the H7 constraint increases 
the sample by a single object; removing both constraints re- 
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Figure 1. Colour versus abs olute magnitude for our current sam- 
ple {blue diamonds) and tlie iPracv et al. I 1I2OO9I ') sample {red tri- 
angles) at redshifts between 0.04 and 0.20. The error bars on the 
Pracy et al. (2009) points give an indication of the uncertainty in 
the k-corrcctions. 



suits in a sample of thirty seven galaxies. The final observed 
targets and their properties are given in Table [1] 

A consequence of selecting an E+A galaxy sample from 
the local volume is that our sample contains intrinsically 
fainter galaxies than those selected at higher redshift. Fig- 
ure[T]shows our sample in a Mr versus q — r color magnitude 
diagram {blue diamonds) along with the lPracv et al.l (|2009l ) 
sample {red triangles). The bright est galaxy in this s ample 
corresponds to the faint -end of the lPracv et al] |2009) sam- 
ple in magnitude and colour. There is a strong correlation in 
the local sample between colour and luminosity in the sense 
that the fainter galaxies are bluer. 



2.2 SDSS data 

The SDSS spectra of our sample are shown in Fig. [2] Three 
galaxies have detectable Ha emission (E+A 2, E+A 3 and 
E+A 5). Two of these galaxies have [Nil] A6583/HaA6563 > 
1 with the other having [NII]A6583/HaA656 3 ~1. Such ra- 
tios a re generally indicative of AGN activity (|Baldwin et al.l 
Il98ll ). although, in this case we have selected galaxies with 
weak Ha emission and underlying stellar populations with 
strong Balmer line absorption. In this circumstance the un- 
derlying absorption can have a large effect on the appar- 
ent ratio of the lines and we show later that these spectra 
are more consistent with a combination of star formation 
and AGN activity. The spectra all have strong higher order 
Balmer lines (as selected) and blue continua - E+A 5 has 
a significantly redder continuum than the remainder of the 
sample but still has an integrated colour that is quite blue 
with g-r ^ 0.78. 

Colour SDSS images of the sample are shown in the 
leftmost column of Fig. [3l The galaxies in the sample gen- 
erally have early type morphologies which is typical of E+A 
samples. There are two exceptions: E+A 1 which has an ir- 
regular morphology and is also the bluest (youngest stellar 
population) galaxy in our sample and E+A 3 which is a late 
type disk galaxy. 
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Figure 2. Spectra of our sample of E+A galaxies from the Sloan Digital Sky Survey data release 7. 
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Table 1. Target galaxies and their properties 



SDSS object ID 


ID 






r 

(mag) 


g ~ V 


z 


Mr 

(mag) 


scale 
(kpc/arcsec) 


H<5 

A 


Morphology and environment notes 


587725073921278177 


E4 


-A 


1 


14.92 


0.20 


0.005 


-16.7 


0.10 


-6.1 


Irregular morphology with lower mass 
starbursting companion 5 kpc away. 


Oc5 1 ( ZOuoOo^ 1 41:01001 




A 


9 


IZ.oO 


U-0c5 


fi ^\(^A 

U.UU'rt 


-10.c5 


u.uo 




IN olOD. isoiateci eariy xype. 
Member of the SAURON sample. 


588017724397912097 


E4 


-A 


3 


15.26 


0.49 


0.009 


-17.7 


0.18 


-5.8 


NGC 3976A. Late type disk located ^50 kpc from 
much more massive NGC 3976. 


588017730298576932 


E4 


-A 


4 


14.50 


0.42 


0.002 


-15.1 


0.04 


-5.4 


In the Virgo cluster. No evidence of interaction. 


588848898841903109 


E4 


-A 


5 


12.80 


0.78 


0.007 


-19.6 


0.14 


-6.0 


NGC 4418. Semi-isolated early type. Starbursting 
companion ~17kpc away. 


587729779053887547 


E4 


-A 


6 


12.98 


0.60 


0.005 


-18.7 


0.10 


-6.3 


Isolated early type 


587726033341513813 


E4 


-A 


7 


14.67 


0.56 


0.007 


-17.7 


0.14 


-6.4 


Isolated early type 



Notes. Basic properties of our sample. Column 1 is the SDSS object ID and Column 2 is our simplified name. The remainder of the 
columns are: SDSS r-band magnitude; g — r colour; redshift; absolute r magnitude; angular to physical scale; H<5 absorption line 
equivalent width from the SDSS integrated spectrum; and a short description of the morphology and local environment from visual 
inspection of the SDSS imaging for each galaxy and the surrounding region. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Observations 

Our data were obtained using th e Wid e Field Spectro- 
graph (WiFeS; lOopita et all l2007l . l20ld ) on the Mount 
Stromlo and Siding Springs Observatory's 2.3-m telescope. 
The WiFeS is an image-slicing integral field spectrograph 
with a 25 arcsecond by 38 arcsecond field-of-view. The spa- 
tial pixel size is 1 arcsecond square (after on-chip binning) 
and the spectra have a resolution of R~3000 and cover 
from ~3500 A to ~9000 A. The WiFeS is designed for Nod 
and Shuffle observations and we used the instrument in this 
mode. This mode of observation would not normally be used 
for observations concerned primarily with the blue end of 
the optical spectrum. However in this case our objects filled 
the entire field-of-view so classical sky subtraction was not 
possible. The nod and shuffle technique was chosen over ob- 
servations of offset sky fields (which would have been more 
efficient) to minimize systematics. The WiFeS CCD is read 
out by four separate amplifiers and as such the four quad- 
rants have different characteristics. We positioned our tar- 
gets away from the centre of the WiFeS field-of view in the 
bottom right quadrant (see Fig. [3| . This was done to avoid 
the galaxy centre being located where the four read-out sec- 
tions meet. The observations were obtained over a 5 night 
dark-time run in March 2010. Table [5] summarizes our ob- 
servations. 



3.2 Data reduction 

The data we re reduced using th e WiFeS IRAF data reduc- 
tion package l|Dopita et al.|[2O10l ) . The raw frames were bias- 
subtracted using the overscan region and trimmed. Since the 
data were taken in nod-and-shuffie mode, sky subtraction is 
performed by shifting the image so that the sky regions are 
moved to the corresponding primary region and subtracted 
from the original image. Flat-fielding was performed utiliz- 
ing quartz-halogen lamp exposures. At this point the indi- 
vidual slice regions were extracted into separate image ex- 



Table 2. Summary of observations. 



ID 


Exposure time (s) 


N 


PA 


seeing 




object sky 


(read outs) 




(arcsec) 


E4-A 1 


2880 2880 


2 


155 


1.7 


E4-A 2 


1200 1200 


1 


50 


1.7 


E4-A 3 


8040 8040 


4 


140 


3.0 


E4-A 4 


6600 6600 


4 


180 


2.0 


E4-A 5 


3600 3600 


3 


60 


1.7 


E4-A 6 


1440 1440 


1 


170 


1.6 


E4-A 7 


3600 3600 


3 


175 


1.7 



Notes. From left to right: Object ID; total exposure times on 
object and on sky; the number of readouts; the position angle of 
the WiFeS IFU; and an estimate of the seeing during the 
observation. 



tensions. The wavelength solution was derived and applied 
using CuAr arc lamp exposures. The slices were then com- 
bined into a single three dimensional data cube. 

3.3 Adaptive binning 

Within each data cube there is a large variation of the spec- 
trum signal-to-noise ratio achieved depending on its spa- 
tial position. For the most part the signal-to-noise ratio 
drops off rapidly moving away from the bright galaxy cen- 
tre. In general, the signal-to-noise ratio of a 1 arcsecond 
spaxel in the central regions is sufficient for quantitative 
analysis (SNR~10 A~^). However, in the outer regions single 
1 arcsecond spaxels do not provide sufficient quality spectra. 
In order to increase the SNR in the outer parts and to ob- 
tain more consistent spectral quality across the field-of-view 
we bin the data spatially using the V oronoi spatial binning 
method of ICappellari fc CopinI (|2003l) . This method, adap- 
tively bins spatial pixels in such a way as to obtain approx- 
imately the same signal-to-noise ratio in each final spatial 
bin. This results in smaller spatial bins in the central region 
and progressively larger bin sizes moving outward. A differ- 
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ent final signal-to-noise ratio per adaptive bin was chosen for 
each galaxy depending on the quality of the data; the final 
signal-to-noise ratios ranged from ~6 to 12 per spaxel. 

3.4 Recession velocities 

We measure the recession velocity of each spatially binned 
spectrum by fitting them with rest frame template spectra. 
The template spectra used were from the si ngle-age single- 
met a llicity stellar population models of IVazdekis et al.l 
( 2007i) which are constructed from empirical stellar spectra 
( Sanchez-Blazguez et al.ll2006l ). The stellar population ages 
used in the models ranged from ~100 Myrs to ~13 Gyrs. The 
spectral resolution of the science data is degraded to match 
that of the templates. The fitting is t hen performed using the 
penali zed pixel fitting algorithm of ICappellari fc Emselleml 
l|2004 ) which fits the spectrum using a combination of input 
template spectra and simultaneously fitting for the recession 
velocity and velocity dispersion. 

3.5 Line indices 

We measure line equivalent widths from each spectrum us- 
ing the flux summing technique. This uses three wavelength 
bands to define the equivalent width; one centred on the 
line itself and two flanking bands on either side of the line 
to define the expected continuum level by linear extrapola- 
tion. For lines at the blue end of the spectrum, including 
H(5, H7 and H/3, we measure line strength equivalent widths 
on the Lick/Intermediate D i spersion Spectrogra ph system 
IWorthev fc Ottavianil 1 19971 : iTrager et all 1 19981 ). Prior to 
measuring the indices the science spectra are convolved with 
a wavelength-dependent Gaussian to the Lick resolution 
(~9 A). For lines red- ward of the Lick indices, in particular 
Hq and [NIIJA6583 that do not have standard definitions, 
we again use the fiux-summing technique with appropriate 
line and continuum bands. 

4 RESULTS 

4.1 Balmer line gradients 

In the third column of Fig. [3] we show the equivalent width 
maps for the H5, H7 and lines combined as a straight 
average. The trends in these maps are also present in each 
of the individual line maps but combining the lines results 
in greater signal-to-noise. The equivalent width maps can 
be compared with the distribution of overall galaxy light 
shown in column 2 by collapsing the light through the orig- 
inal unbinned IFU elements along the wavelength direction 
and column 1 which shows an SDSS colour image with the 
IFU field-of-view superimposed. In every case except E+A 
1 there is strong enhancement in the Balmer line absorption 
equivalent width concentrated in the galaxy core on scales 
< 1 kpc. E+A 1 has a uniform very strong H(5 signature 
spread across the entire extent of the galaxy. It also differs 
to the rest of the sample in other respects: it is the only 
irregular system and easily the bluest galaxy in the sample. 
As mentioned in Table[T]E-|-A 2 (NGC 3156) is a member of 
the Spectroscopic Areal Unit For Research on Optical Neb- 
ulae (SAURON) sample (|de Zeeuw et all 120021 ') enabling us 



to com pare our Balmer absorp tion line map with the II/3 
map of iKuntschner et alj and we note that the maps 

are in good agreement. In the last column of Fig. Owe have 
azimuthally binned the data before measuring line indices 
to construct radial Balmer line gradients. The Balmer line 
enhancement seen in the two dimensional maps are seen as 
negative Balmer line radial gradients. These radial profiles 
of the Balmer line equivalent width are similar to th o se pre - 
dicted from merger models c.f. Fig. 9 of iPracv" et all (|2005h . 
The observed gradients are confined to a smaller radius than 
the models (which show enhancements with a negative gra- 
dient over ~2 kpc) but this is expected given the models are 
for signiflcantly brighter galaxies which should have a larger 
central starburst region. 

In Figure 2] we show an age-metallicity diagnostic dia- 
gram for our galaxies using the 115 and C24668 lines. The 
line ratios are shown for the radial data with the smallest 
red diamond representing the innermost radial point and in- 
creasing radii represented by increasing larger symbol sizes. 
The age-metallicity grids are from the single stellar popula- 
tion models of Thomas ct al. (2003, 2004) and assume solar 
abundance ratios. The Balmer line radial gradients can be 
seen as luminosity weighted age gradients in this plot. In 
general, the smaller inner radial data points have younger 
ages than those at larger radii. The inner radial regions have 
s^ges < 1 Gyr consistent with what is expect ed for E-l-A 
galaxi es and similar to the global ages of the iPracv et al.l 
|2Q0i) E-l-A sample which are overlaid as green squares. 

4.2 Mass fractions 

We estimated the mass in recently formed stars from the 
best fltting stellar population templates as described in Sec- 
tion 3.4. In this case we constructed a spatially integrated 
spectrum for each of our sample by summing over all spax- 
els within a circular aperture of radius 11 arcseconds. We 
fitted the spectra with a combination of a single old stellar 
population of age 13 Gyr and a set of young stellar pop- 
ulation models with ages ranging from 100 Myr to ~ 1 Gyr 
(all model stellar populations assumed solar metallicity) and 
calculated the mass fraction of the young stellar population. 
This results in burst mass fraction estimates of between 6 
and 12 per cent for all galaxies in our sample except E-f-A 1 
which is best fitted with 100 per cent young mass fraction i.e. 
we are unable to quantify the old stellar population in this 
galaxy. These mass fractions correspond to the young stel- 
lar population contributing between 51 per cent and 63 per 
cent of the light between 4500A and 5500A. The calculated 
burst mass fractions are not sensitive to the precise age cho- 
sen for the old stellar population since the spectra of old 
stellar populations do not change rapidly with age. For the 
young population there are degeneracies between age and 
burst fraction in the sense that decreasing the age of the 
stellar population has a similar effect as increasing the burst 
fraction. To illustrate this we fitted again allowing the com- 
bination of only two models. First we fitted an old ^ 13 Gyr 
model and a young stellar population of age 890 Myrs which 
resulted in the mass fraction estimates increasing to between 
11 per cent and 17 per cent (E-|-A 1 again returns a 100 per 
cent young mass fraction and is poorly fitted by this combi- 
nation). We also fitted the combination of a ~ 13 Gyr model 
and a younger 500 Myr template which resulted in the range 
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Figure 3. One object per row - from left to right: SDSS image with the size, orientation and position of the WiFeS field-of-view 
superimposed; reconstructed image by collapsing the WiFeS data in the spectral direction; Balmer line equivalent width maps for 
H5 + H7 + H/3. The overlaid contours are the continuum light profiles; radial Balmer line gradients based on azimuthal binning of the 
WiFeS spectra. Note that the horizontal axis scale in the final column changes. 
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Figure 4. Age metallicity diagnostic diagrams using the H<5 and C24668 lines. The radial annular equivalent width measurements are 
displayed as red d iamonds with symbol sizes increasing with galacto-centric radius. The green squares are the global values from the 
iPracv et al"] | |2009|) sample. 



of mass fractions decreasing to ~3-4per cent (E+A 1 has 
a young mass fraction of ~25 per cent under these assump- 
tions). These lower burst fractions by mass still contribute 
significantly to the overall light (between 40 and 50 per cent) 
since the younger stellar populations have a smaller mass-to- 
light ratio. Reducing the age of the young fraction reduces 



the burst mass fractions further but results in progressively 
poorer fits (again the exception is E+A 1 which is better fit- 
ted by a younger population). Overall, the recently formed 
stars account for of order 10 per cent of the total stellar mass 
but contribute more than half the total light at ~5Q0QA. 
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Figure 5. Shows the effect of changing the redshift each galaxy 
is observed at on the H(5 equivalent width measured by an SDSS 
fibre. We plot H(5 equivalent width measured in different apertures 
versus the redshift that galaxy would be at to make the aperture 
3 arcseconds - the diameter of the SDSS fibres. 



4.3 Aperture effects 

The presence of such strong radial gradients in the Balmer 
hnes means the classification of these objects are strongly de- 
pendent on aperture effects. The measured equivalent width 
of the Balmer lines will depend on the physical scale probed 
which in turn depends on both the distance of the object 
and the angular size subtended by the fibre. At a given dis- 
tance a smaller fibre will measure a larger absorption line 
strength since the light is dominated by the central region 
with a large equivalent width, as the fibre size is increased 
the light becomes more contaminated by the outer regions 
where the absorption equivalent width is small. Likewise, for 
a fixed aperture the Balmer line equivalent width measured 
will decrease with distance as the angular size of the galaxy 
decreases and more of the outer regions of the galaxy are 
within the fibre aperture. This means, independent of ap- 
parent magnitude constraints, galaxies that fit the equiva- 
lent width selection criteria at low redshift may not if they 
were observed at higher redshift. This is illustrated for our 
sample in Fig. [5] where we have measured the H5 equiva- 
lent width for different sized apertures. The horizontal axis 
shows the effective redshift the object would be at to sub- 
tend a three arcsecond aperture (which is the fibre aperture 
of the SDSS) i.e. we have measured the H5 equivalent width 
in our IFU data with an aperture that is the same physical 
size as a 3 arcsecond aperture would be at that redshift. This 
demonstrates how the H5 equivalent width measured would 
decrease for our objects if they were observed at higher red- 
shift. The horizontal solid black line indicates the US equiv- 
alent width selection limit for our sample. All of our sample 
(except E+A 1 which does not display a radial gradient) 
would have H(5 equivalent widths below the selection limit 
if observed at z > 0.06. 



4.4 Emission lines 

As mentioned in Section 2.2, three of the galaxies have 
integrated spectra that have Ha in emission. This emis- 
sion is weak, by selection, since we required EW(Hq)>- 
2.5A. In ah three cases the ratio [Nil] A6583/HaA6563 >1 
which is indicative of AGN activity rather than star forma- 
tion fe.g [Baldwin et al. l ll98ll : IVeilleux fc Osterbrockl [l987l : 
iKewlev et al.l 2001 ). However, as a result of our selection 
of spectra with both weak or non-existent Ha emission 
and a stellar continuum with strong Balmer lines in ab- 
sorption, we therefore expect strong Ha absorption in the 
stellar continuum underlying the emission line. Hence the 
[Nil] A6583/HaA6563 ratio is significantly overestimated. 

In Fig. [6] we show spatially resolved emission line maps 
for the [Nil] and Ha lines in columns 1 and 2, respectively. 
The three galaxies that show Ha in emission in the inte- 
grated SDSS spectra (E+A 2, E+A 3 and E+A 5) in Fig. 
[2I have their emission confined to the central region. In the 
case of E+A 3, there is convincing emission in only a single 
spaxel in the galaxy centre and most of the galaxy has Ha 
in absorption. The rest of the sample shows no evidence for 
emission at larger galactic-centric radii that would not have 
been contained in the 3 arcsecond SDSS aperture. The weak 
values of Ha emission in E+A 1 are essentially due to noise, 
and inspection of the spectra (which are poor quality for this 
galaxy) shows little evidence for emission or absorption. 

In column 3 of Fig. [6] we attempt to correct for the issue 
of stellar absorption on the Ha emission line index measure- 
ments. This is done by measuring the Ha line equivalent 
width from the spectrum of the best fitting stellar popu- 
lation model to the blue arm data. This gives an estimate 
of the strength of the underlying Ha absorption due to in- 
tegrated light from the stellar atmospheres alone. This can 
then be subtracted from the observed line equivalent width 
to obtain an estimate of equivalent width of the emission. 
Maps with this correction to the Ha emission are shown in 
column 3 of Fig. [G] The central Ha emission is still clear in 
E+A 2, E+A 3 and E+A 5 with stronger absolute equiva- 
lent widths. In addition E+A 4 and E+A 6 also show a weak 
enhancement of emission in the galaxy centre and E+A 7 
shows greater emission in the galaxy outskirts. E+A 1 has 
uniform emission after correction of ~5A, which is the result 
of the observed spectra showing little emission or absorp- 
tion at Ha but strong absorption in the bluer Balmer lines 
causing a large correction. The outer regions of the rest of 
the sample show ~1 to 2 A of emission line filling indicat- 
ing there may be small amounts of ongoing star formation 
or that the fit to the higher order Balmer lines is slightly 
overestimating the Ha equivalent width. 

In the final column of Fig. [H] we plot the ratio of 
[NII]A6583/HaA6563 where we have used the corrected 
Ha value. We only plot spaxels which have EW(Ha) and 
EW([NII]A6583) > O.2A and a signal-to-noise ratio in both 
these lines greater than 2. The most common diagnostic 
plot for classifying emission line galaxies is the BPT dia- 
gram (Baldwin et al. 1981) and uses the position of a galaxy 
in the two dimensional [0III]/H/3 versus [Nil] /Ha plane. 
However, this diagnostic cannot be used for many weak line 
galaxies or low signal-to- noise spectra where the H/3 and/or 
[OIII] lines cannot be satisfac torily m easured - which i s 
the case here. For this reason ICid Fernandes et al.l (|201ll ) 
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Figure 6. Emission line strength and diagnostic maps. The first column is the [NII]A6583 equivalent width map; the second column is the 
Ha emission equivalent width map; the third column is the Ha equivalent width map after correction for underlying stellar absorption; 
and the final column shows the [Nn]A6583/Ha ratio (corrected) which can be used to classify the source of the ionizing radiation. The 
contours overlaid are the galaxy integrated light isophotes. 
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developed several more inclusive diagnostics for emission 
line taxonomy, the simplest of which uses the EW(Ha) ver- 
sus [NII]/Hq plane, which is available here. Optimal trans- 
positions of the common division lines for classifying ob- 
jects as star-fo rming, Seyfert or L I NER in the traditional 



BPT diagram jKewley et al.l I2001I: 



Stasihska et al. 



2006; Kcw lev et al 



Kauffmann et all [200^ : 
I2OO6!) are derived by 
Cid Fernandes et a l. (2011). In this plane the standard divi- 
sion b etween star formation and AGN of iKauffmann et al.l 
l|2003tl translates to [NII]A6583/HqA6563 ~0.5 with values 
lower than this being in the star formation regime. The 
hard u pper limit for star forming galaxies of iKewlev et al.l 
l|200l[) translates to [NII]A6583/HaA6563 ~0.8 with galax- 
ies between these t wo values being com posite systems and 
galaxie s above the iKewlev et al.l (|200ll ) value being pure 
AGN (|Cid Fernandes et al.l I2OI1I ). Comparing these val- 
ues with the maps in the final column of Fig. 6 for the 
galaxies with Ha emission shows E-l-A 2 and E-l-A 5 have 
[NII]A6583/HaA6563 ~0.5-0.7 corresponding to the tran- 
sition region whilst E-l-A 3, the only late type disk sys- 
tem in the sample, has a line ratio consistent with pure 
star-formation. A bright point-like continuum source was 
detected in the ce ntre of E-l-A 5 in the FIRST survey 
l|Becker et al.lll995l ) consistent with the presence of an AGN 
in this galaxy. 

Since our spectra extend blue- ward of the [OII]A3727 
line we can check whether the three galaxies with Ha emis- 
sion would have been selected in the sample if the [OII]A3727 
line was available and the normal [OII]A3727 > — 2.5A se- 
lection cut was applied. We measured the [OIIJA3727 line 
equivalent width in a 3arcsecond aperture centred on the 
galaxy to match the SDSS fibre size. The [OIIJA3727 equiv- 
alent width of E-l-A 2 is ~ — 1.9±0.7A satisfying the nor- 
mal selection criterion. E-f A 3 and E-f A 5 have [OII]A3727 
equivalent widths of ~ —2.7 ± 2.3A and —5.2 ± I.3A, re- 
spectively, both slightly above the threshold, although in 
the case of E-|-A 2 it is well within the error. The rest of the 
sample have [OII]A3727 > —2.5 A. and consistent with zero 
at la. 



4.5 Kinematics 

Two dimensional streaming velocity maps are shown in Fig. 
[7] and all the galaxies except E-l-A 1 and E-l-A 3 show ro- 
tation at some level. In the case of E-l-A 3 this is surpris- 
ing since it is a disk galaxy viewed somewhat edge on and 
rotation is to be expected. E-l-A 3 is the faintest and high- 
est redshift galaxy in our sample and was observed in very 
poor seeing and it maybe the case that the data quality 
is insufficient to detect the velocity field. Using a parame- 
ter, Aij, derived from two-dimensional spectroscopy, early- 
type galaxies can be separated into two distinct kinemati c 
categories: fast and slow rotators (|Emsellem et al.l |2007^ . 
The \r parameter involves a luminosity weighted average 
over the kinematic field ob tained from IFU spectroscopy. 
lEmsellem et all (|2007l .[ 201 ll ) found that approximately three 
quarters of the SAURON sample were fast rotators. Like- 
wise, in a sample of eight E-l-A galaxies selected from the 
2dFGRS and followed up with int egral field spectro scopy 
seven were found to be fast rotators l|Pracv et al.ll2009h . The 
prevalence within the E-l-A population of fast rotators ar- 
gues that unequal mass mergers rather than major mergers 



were the dominant progenitors to the E-l-A galaxies since the 
probability of rota ting remnants increas es with increasing 
merger mass ratio ( jBournaud et al. | |2008l ). Altho ugh major 
merg ers can still produce fast rotating remnants (|Bois et al.l 
I2OI 



I 



In Fig. [S] we plot the values of \r for the early type 
galaxies in our sample (not including E-l-A 1 and E-l-A 3 
which are irregular and late type disk systems) versus the 
radius over which it is measured (top panel) and against el- 
lipticity at one effective radius (bottom panel). The stream- 
ing velocities and velocity dispersions used to calculate A_r 
were derived from the spectral fitting described in Section 
3.4. We overlay in Fig.[S]the data from the SAURON sample 
for fast rotators {blue lines and symbols) and slow rotators 
{red Imes and symbols). In all five cases our data overlap the 
fast rota tor regions of thes e plots consistent with previous 
findings l|Pracv et al.ll2009l ). The Ah values for NGC 3156 
(E-f A 2) from the SAURON sample are shown as a thick 
line and large filled diamond in the top and bottom panels 
respectively. These can be directly compared with the value 
for E-f A 2 from our measurements; the derived values of A_r 
are similar with the SAURON study value ^0.1 larger. 

There is a minimum accurately measurable value for 
\r l|Emsellem et al.ll2007l ) which depends on the quality of 
the data. Since Xr depends on the square of the velocity of 
each spaxel, even for data with zero true rotation, random 
noise will induce a non-zero positive value of \r. To deter- 
mine if our derived values of \r are reliable we simulate this 
effect by setting the velocity values in each spaxel to zero 
and adding a random velocity drawn from a normal distri- 
bution with a width determined by the velocity error for 
that spaxel. We use 10000 such simulations for each galaxy 
to estimate the minimum value of \r that can be reliably 
measured. For E-f A 2 and E-f A 5 this value is <15 per cent 
of the measured value and these galaxies can be confidently 
classified as fast rotators. For E-f A 7 this value is ~40 per 
cent and in the cases of E-f A 4 and E-f A 6 the minimum 
measurable values are of order 60 per cent of the measured 
values. In a similar way we calculated the errors on the A_r 
values by randomly varying the streaming velocity and ve- 
locity dispersion of each spaxel and using the scatter in the 
resulting values of \r as an estimate of the error. 



5 DISCUSSION 

By selecting a sample of E-f A galaxies that are considerably 
closer, but less massive, than previously studied systems, 
we have been able to obtain spatially resolved spectroscopy 
of the core (Ikpc) region of E-f A galaxies for the first 
time. Models that form E-f A galaxies through tidal inter- 
actions or mergers suggest that starbursts occur in the core 



core (iBarnes & HernauistI 19961: 


Mihos & HernauistI Il994l: 


iHoDkins et 


all 20091: Snvder et al. 


I2OIII). This region should 


be compact 
[Bckki ct al. 


on scales of ~1 kpc, ( 


Barnes & HernQuistlll996l: 


2005; HoDkins et all 


2009|) and contain strong 



stellar population g r adien t s observable as Ba l mer line gra- 
dients l|Pracv et al.l l2005l : iBekki et al.l l2005l : ISnvder et all 
l201lD . The primary motivation of our observations was to 
test these predictions by acquiring spatially resolved spec- 
troscopy that resolved this region. 
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Figure 7. Shows the stellar streaming velocity maps for our sample. E+A 1 and E+A 3 have no detectable rotation and E+A 4 has 
only marginal evidence for rotation. The remaining galaxies all show clear rotation fields. The black contours overlaid are the galaxy 
integrated light profiles. 



Previous stu dies, that used colours from imaging 
jYang et al. 2008") or spectroscopy that docs not resolve the 
core l|Norton et al. 2 001 ; Chiling arian ct ah 2009), have sug- 
gested that the post-starburst region is centrally concen- 
trated but not necessarily confined to the galaxy core. Other 
groups have been unable to detect any evidence for a cen- 
tral concentration in their observations (|Yagi fc Gotoll2006l : 
iPracv et al]|2009h . The most distant object in our sample 
has a redshift of z ~0.009 corresponding to a physical to an- 
gular scale of ~0.18 kpc per arcsecond which allowed us (for 
this worst case object) to easily spatially resolve the cen- 
tral 1 kpc region even in moderate seeing conditions. The 
resulting Balmer line equivalent width maps confirm that 
the post starburst region is centrally concentrated and con- 
tained within <lkpc (see third column of Fig.|3} consistent 
with model expectations. The radial profiles of the Balmer 
line absorption maps which reveal a steep negative Balmer 
line gradient are consistent with the model predictions for 
E+A galaxies pro duced from mergers or tidal interactions 
IjPracv et al. 1 12005!). The exception is E+A 1 which has uni- 
formly strong Balmer line absorption and no radial gradient. 
E+A 1 differs from the rest of the sample in other respects: 
it is the only galaxy with an irregular morphology, it is the 
bluest galaxy in the sample by a considerable margin and 
has the strongest Balmer line equivalent width values imply- 
ing a younger post-starburst age. After emission line filling 
is taken into account it has stronger residual Ha emission 
than the rest of the sample indicating more ongoing star- 
formation which is distributed over the entire galaxy. E+A 
1 may represent an earlier evolutionary stage than the other 
galaxies in the sample or a different class of object alto- 
gether. 



It has been pointed out bv lSnvder eraD \2Ql± that the 
presence of strong radial Balmer gradients in E+A galaxies 
means that the galaxies selected by a given equivalent width 
criterion depend on the galaxy distance since a fixed diam- 
eter fibre will subtend a greater physical distance at higher 
redshift. For centrally-concentrated Balmer gradients this 
implies that for a given selection criterion a more extreme 
population will be selected at higher redshift. This is due to 



the aperture at higher redshifts including more light from 
larger galacto-centric radii where the Balmer line absorp- 
tion is weaker. At the median redshift of our sample, the 
3 arcsecond fibre of the SDSS subtends just ~300pc whereas 
at the median redshift of the SDSS (z ~0.1) a 3 arcsecond 
fibre subtends ~5kpc. For our galaxies the strong radial 
gradients and strong 'E+A levels' of Balmer line absorp- 
tion are contained within the central few hundred parsecs. 
Similar galaxies at higher redshift would not necessarily be 
classified as E+A galaxies. Selecting galaxies based on mea- 
surements from a central fibre also biases us against selecting 
galaxies which have positive H5 gradients (central deficits) 
which are expecte d from the sudden truncation of star for- 
mation in a disk (jPracv et al.ll2005l ). Our sample contains 
on average less luminous galaxies and also objects with less 
extreme overall absorption line strengths than many of the 
previous E+A samples in the literature. 

The morphologies of the galaxies in our sample, like 
E+A samples in general, are mostly early types with the ex- 
ceptions being one late-type irregular galaxy and one late- 
type disk galaxy. Most of the early type galaxies are isolated 
systems with no interacting neighbor present - see Table [1] 
for a list of morphologies and environments for the sample. 
This isolation rules out ongoing interactions as the cause 
of the post-starburst signature and argues in favor of these 
galaxies being merger remnants where the progenitors have 
already coalesced and relaxed. The kinematics of the early 
types in our s ample are all consistent with being fast rotators 
as defined bv lEmsellem et al] l|2007l ). This suggests that the 
mergers producing the E+A galaxy are often unequal mass 
or minor mergers since mergers with mass ratios around 
unity will more comm only result in a non-rotating remnant 
(|Bournaud et al.l l 2008|), although, th is does not strictly rule 
out major mergers (Bois et al.' 2011). This same conclusion 
was reached by Pracy ct al. (2009) for E+A galaxies selected 
from t he 2dFGRS and is in a greement with the recent find- 
ings of ICrockett et all (|201ll ). The two late type galaxies in 
our sample (E+A 1 and E+A 3) both have near neigbours 
which are very close in the case of E+A 1(^5 kpc) or much 
more massive in the case of E+A 3 (see Table [1] for details) 
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6 CONCLUSIONS 

We have selected a new local sample of seven E+A galax- 
ies from the SDSS with a redshift limit of 2: = 0.01 and 
obtained integral field spectroscopy with the WiFeS instru- 
ment on the ANU 2.3-m telescope. The low redshift of our 
sample allowed for much better physical scale resolution in 
the core region of the E-f-As than previous studies. This 
allowed us to unequivocally verify that the post starburst 
regions are centrally concentrated, consistent with expecta- 
tions. In summary our main conclusions are: 

• Six of the seven E-l-A galaxies and all galaxies with regu- 
lar morphologies have a centrally-concentrated young stellar 
population and negative Balmer line gradients on scales of 
<lkpc. The exception is the only irregular galaxy in the 
sample which shows uniform Balmer absorption strength 
and stronger overall absorption compared to the other galax- 
ies. 

• The sample is dominated by isolated early type galaxies. 
All the early type galaxies show some level of rotation and 
are consistent with the fast rotator population. 

• In the early type galaxies the combination of a centralized 
young population and isolation argues in favor of a merger 
origin where the progenitors have already coalesced. 

• The two late type galaxies in the sample have nearby com- 
panions and their spectral signature could be explained by 
tidal interaction. 



Figure 8. Top panel: the parameter versus the radius over 
which it was measured for the early type galaxies in our sample. 
We measured Xft over a circular aperture which is 85% covered 
by the IFU data and correct the radius to an effective circular 
radius which has the s ame areal coverage. Th e tracks for from 
the SAURON sample jEmsellem et al.ll2007l) are overlaid as blue 
dashed lines for the fast rotators and as red solid lines for the 
slow rotators. Bottom panel: the A^ parameter versus ellipticity 
measured at one effective radius. Again the SAURON data arc 
overlaid for fast {blue diamonds) and slow {red squares) rotators. 
The SAURON data for NGC 3156 (E+A 2) are plotted as a thick 
line (top panel) and filled large diamond (bottom panel) 
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and for these two galaxies the E+A signature could be ex- 
plained by tidal interactions rather than merging. 

There are detectable emission lines in two of the early 
type galaxies as well as in the late type disk galaxy. The 
emission is centrally concentrated and absorption-corrected 
line ratios imply that in the early type galaxies this emission 
is, at least in part, due to AGN activity. The line ratios in 
the disk galaxy are more consistent with pure star forma- 
tion. There are o ther known examples of low power AGN 
in E+A galaxies (|Liu et al.|[2007l : IChilingarian et al.|[2009l ) 
and it has been proposed that nuclear activity, in addition 
to supernova feedback, contributes to the truncation of star- 
formation in E+As tKavirai et al., ,2007i ). However, simula- 
tions have shown that the eff ect on the post-starburst phase 
from AGN feedback is weak (|Wild et al.ll2009l : ISnvder et all 
I2OIII ). with the dominant effect being an increase in the ob- 
served post-starburst signature as a result of the expulsion 
of obscuring dust. 
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